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ABSTRACT 

Wide area X-ray and far infrared surveys are a fundamental tool to investigate the link between AGN growth and 
star formation, especially in the low-redshift universe {z < 1). The Herschel Terahertz Large Area survey (H-ATLAS) 
has covered 550 deg^ in five far-infrared and sub-mm bands, 16 deg^ of which have been presented in the Science 
Demonstration Phase (SDP) catalogue. Here we introduce the XMM-Newton observations in H-ATLAS SDP area, 
covering 7.1 deg^ with flux limits of 2 x 10“^^, 6 x 10“^^ and 9 x 10“^^ erg cm“^ in the 0.5-2, 0.5-8 and 2-8 
keV bands, respectively. We present the source detection and the catalogue, which includes 1700, 1582 and 814 sources 
detected by EMLDetect in the 0.5-8, 0.5-2 and 2-8 keV bands, respectively; the number of unique sources is 1816. We 
extract spectra and derive fluxes from power-law fits for 398 sources with more than 40 counts in the 0.5-8 keV band. 
We compare the best-fit fluxes with the catalogue ones, obtained by assuming a common photon index of F = 1.7; we 
find no bulk difference between the fluxes, and a moderate dispersion of s = 0.33 dex. Using wherever possible the 
fluxes from the spectral fits, we derive the 2-10 keV Log A^-Log S, which is consistent with a Euclidean distribution. 
Finally, we release computer code for the tools developed for this project. 
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1. Introduction 


Over more than one decade, there has been growing ev¬ 
idence for a coeval growth of galaxies and their cen tral 
black holes (see review bv [Alexander fc HickoxI r2012f ). A 
tight correlation between the masses o f the black hole and 
the galaxy bulge has been found (e.g. Ferrarese fc Merritt 
2QQQI : iGebhardt et al.l l2QQQl: IZubovas fc Kind l2Q12f ) . Theo¬ 
retical models suggest that feedback processes are at wor k 
to set the link (|Di Matteo et al.ll2QQ5l:lHQDkins et al.ll2QQ6[ ): 
some models also s uggest a fundament al role of mergers 
in setting up AGN ([Hopkins et al.l[200^ . Nuclear obscura¬ 
tion and intense star forma tion may characterize the ini¬ 
tial phases of AGN activity ([Silk fc Reeslll99^: iMenci et al.l 
2 Qq 4 iLamastra et al.ll2QT^ . 


AGN growth seems to happen in two major modes: 
the radiative and the kinetic mode, the first operating 
close to the Eddington limit and with high radia tion ef¬ 
ficien cy, the latter at lower rates (see review by iFabianI 
l2Q12f ). Similarly, galaxies build their stellar mass either 
through starburst episodes (with star formation happen¬ 
ing on short timescales), or through secular star formation. 
The growth of both AGN and galaxy should ultimately be 


* Based on observations obtained with XMM-Newton, an ESA 
science mission with instruments and contributions directly 
funded by ESA Member States and NASA. 


driven by the supply of col d gas ([Kanffmann fc HeckmanI 
I2nnfll: iMiillanev et, ai jl2nT^ . 

The fraction of galaxies hosting an AGN increases 
with far-infrared l uminosity or s t ar formation rate (SFR ) 
(|Kim et al.l 1 19981: IVeillenx et al.l 1 19991: iTran et ^ l2QQl[ ). 
reaching the 50-80% among Luminous InfraRed Galaxies 
LIRG) and UltraLum i nous InfraRed Galaxies (ULIRG) 
Alexander et al.l 20081: Lehmer et al.l |201(]I: iNardini et alJ 


20101: iNardini fc Risaliti 1^111: lAlonso-Herrero et al.ll2012l: 

Ruiz et al.ll2013f ): IR and X-ray observations are the key 


for the identification of AGN. 

It has been suggested that the specific SFR (i.e., the 
SFR divided by the stellar mass of the galaxy) may 
also be involved in the AGN growth/star formation link, 
because a correlation was observed betwe en AGN lumi¬ 
nosity and specific SFR ([Lutz et al.l l2010f) : however the 
latter correlation seems to hold only for the most ac¬ 
tive systems at re dshift z > 1 ([Mnllanev et al.l l2012bl: 
iRovilos et al.|[2QT^ . Similar uncertainties shroud a possible 
corre lation between AGN lum i nosity and nuclear obscura¬ 
tion ( ^orgakakis et al.| [2006l: iRovilos fc GeorgantoponlosI 
l2QQ7l: Trichas et al.ll2QQ9l: IRovilos et aOEo^l^ ) and the effec¬ 
tiveness of colour-magnitude diagrams to insp ect the evo¬ 
lutionary status of the host galaxies of AGN (iBrnsa et al.l 


_ t ga _ 

l2QQ9l: ICardamone et al]l2QlQl: [Pierce et ^l2QlQf) . 

The situation in the local universe is mostly unclear, 
as the aforementioned studies were based on deep, pencil- 
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Table 1. XMM-ATLAS observations. The columns show: observation ID; date; total exposure times (ks) for the MOSl, MOS2 
and PN cameras; exposure times (ks) after high background filtering for the MOSl, MOS2 and PN cameras; astrometry corrections 
applied (arcsec). 


Obsid 

Date 

MOSl 

MOS2 

PN 

MOSl clean 

MOS2 clean 

PN clean 

ARA 

ADec 

0725290101 

2013-05-05 

113 

113 

no 

93 

96 

71 

— 

— 

0725300101 

2013-05-07 

113 

113 

no 

101 

100 

61 

-0.44 

-0.47 

0725310101 

2013-05-21 

no 

no 

112 

99 

100 

94 

-1.43 

+0.41 


beam surveys. Wide area surveys are hence needed to probe 
larger volumes of the low-redshift universe, and build siz¬ 
able samples of rarer objects. 

The Herschel Terahertz Large Area survey (H-ATLAS) 
is the largest Open Time Key Project carried ou t with the 
Herschel Space Observatory (|Eales et al.l l2QlQh . covering 
^ 550 deg^ with both the SPIRE and PACS instruments 
in five far-infrared and sub-mm bands (100, 160, 250, 350 
and 500 jam). The 250 jam source catalogue o f the Science 
Demo nstration Phase (SDP) is presented in iRigbv et al.l 
(1201 It) . covering a contiguous area of ^ 16 deg^ which lies 
within one of the regions ob served by the Galaxy And Mass 
Assembly (GAMA) survey (|Driver et al.ll2009l:lBaldrv et al.l 
1201011 . 

In this paper, we present XMM-ATLAS i.e. the XMM- 
Newton observations and source catalogue in the H-ATLAS 
area. XMM-Newton observed 7.1 deg^ within the H-ATLAS 
SDP area, making the XMM-ATLAS one of the largest con¬ 
tiguous areas of the sky with both XMM-Newton and Her¬ 
schel coverage. 

The only other wide X-ray survey with size and Herschel 
cover age comparable to XMM-ATLAS is the 11 deg^ XMM- 
LSS ([Pierre et al.ll2003: IChiaopetti et al] 120131 1 which was 
observed with SPI RE and PACS by the HerMES project 
([Oliver et al.ll2012[ l. Th e Stripe-82 survey cov ered 10.5 deg^ 
with XMM-Newtor^ ([LaMassa et al.l 120131 1. though the 
pointi ngs are not contig uous and it only has SPIRE cov¬ 
erage ([Viero et al.ll2014l l. The 2 deg^ COSMOS survey , ob¬ 
served by both XMM-Newton (ICappelluti et aLll^OPOl l and 
Chandra ([Elvis et al.ll2009l: [C ivan dl2ni3h. was als o covered 
by HerMES with SPIRE only ([Oliver et all 1201 2[ 1. 

Eor the bright sources in the XMM-ATLAS catalogue, 
we extract spectra and derive fluxes from the spectral fits. 
The latter fluxes are used to build the Log A^-Log S', thus 
removing as much as possible any bias which might come 
from using a single count-rate-to-flux conversion factor as 
done for the catalogue. 

In Sect. [21 we describe the observations and data re¬ 
duction; in Sect. |3l we illustrate the source detection; in 
Sect. lUwe present the source catalogue in the 0.5-2, 0.5- 
8 and 2-8 keV bands; in Sect. [5] we derive spectra for the 
bright sources and compare the fluxes from the spectral fits 
with those from the catalogue; in Sect, ^we compute the 
Log N-Log S using, where available, the fluxes from the 
spectral fits. Einally, in Sect. [7| we present our conclusions. 

2. Observations and data reduction 

The XMM-ATLAS field is centred at 9h 4m 30.0s 
+0d 34m Os, and covers 7.101 deg^, with a total exposure 
time of 336 ks (in the MOSl camera). The observations 
were performed in mosaic mode, i.e. shifting the pointing 

^ Reaching 16.5 deg^ when considering both XMM-Newton and 
Chandra data. 



Fig. 1. Distribution of the exposure times in the XMM-ATLAS 
mosaic (bin size: 1 ks). The times are the average between 
MOSl, MOS2 and PN, for the 0.5-8 keV band; the vertical axis 
shows the fraction of pixels with a given exposure. 

by 15' every ^ 10 ks. A total of 93 pointings were done, 
divided in 3 obsids of 31 pointings each. 

The SAS (version 13.0) tools emproc and epproc were 
used to produce a single event file per obsid per camera. 
Such an event file can be directly used to obtain images 
and exposure maps of the part of the mosaic covered in 
the obsid. However, the file needs to splitted in order to 
produce images and exposure maps of individual pointings 
(see Sect. (Ml). 

We extracted light curves in the 10-13 keV interval with 
a 100 s bin size to check for high-background periods, which 
were identified and re moved by doing a 3(j -clipping of the 
light curve, as done in iRanalli et aLl ((20131 . hereafter R13) 
for the XMM-CDES. After cleaning, the total exposure is 
293 ks (MOSl). 

The obsid, dates and exposure times before and after the 
high-background cleaning, are listed in Table [U note that 
these times refer to the mosaic-mode obsid. The average 
exposure time for any location in the final mosaic is 3.0 
ks (in the 0.5-8 keV band and after the high-background 
cleaning); the maximum exposure is 11 ks. A histogram of 
the exposure times is shown in Fig. [H 

To check if astrometric corrections were needed, we 
performed an initial source detection run with the SAS 
ewavelet program (with a 5cr threshold) on the PN data 
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Fig. 2. Left: image of the entire mosaic of the XMM-ATLAS observations in the 0.5-8 keV band. The greyscale wedge shows the 
photon counts (truncated; the actual maximum is 117 counts). A Gaussian smoothing with 20^^ radius has been applied to enhance 
the sources. Some areas with longer exposure times and/or higher background are visible. Right: exposure map in the 0.5-8 keV 
band. The greyscale wedge (maximum value: 33.9 ks) shows the sum of the MOSl, MOS2 and PN exposure times in s, corrected 
for vignetting and bad pixels. 


of each of the three obsids, and cross-correlated the result¬ 
ing lists with the sampl e of QSO from the S loan Digital 
Sky Survey (SDSS DR7: IS^neider et al.ll2QlQl) . The search 
radius was 15", and most matches were found within 5". 
We found that obsid 0725290101 (for which 26 matches 
were found) needed no correction; obsid 0725300101 (35 
matches) had a shift of 0.58" and obsid 0725310101 (39 
matches) of 1.14". Event and attitude files were corrected 
(see Tahiti]) before proceeding further. 

Images for the MOSl, MOS2 and PN cameras were ac¬ 
cumulated for each obsid and for each of the 0.5-2, 2-8 
and 0.5-8 keV bands with a pixel size of 4" and summed 
together. Some energy int ervals corresponding to k nown in¬ 
strumental spectral lines (|Knntz fc Snowdei]120081 see also 
R13) were excluded: 1.39-1.55 keV (Al); 1.69-1.80 keV (Si, 
MOS only); 7.35-7.60 and 7.84-8.28 keV (Cu comple:^^, PN 
only). Images of the entire mosaic were obtained by sum¬ 
ming the images and exposure maps from the three obsids; 
those relative to the 0.5-8 keV band are shown in Fig. [21 

3. Source detection 

Sources were detected with a two-stage process, with a first 
pass at low significance to get a list of candidate detec- 

^ The same complex also includes a line in the 8.54-9.00 keV 
interval, which is outside of the bands considered here. 


tions, and a second pass to rise the significance threshold 
and derive accurate source parameters. Between the two 
passes, and because the second pass needs an input cata¬ 
logue, we identify the sources detected in more than one 
band. The detec t ion m ethod is based on the one used by 
iCaPDellnti et al.l (|2009[ ) and R13, which is a variant on the 
standard XMM-Newton detection procedure, adapted to 
the much larger area of XMM-ATLAS. 


3.1. First detection pass 

In the first pass, the SAS wavelet detection program 
e wave let was run separately on each of the 0.5-2, 2-8 and 
0.5-8 keV images of the entire mosaic, with a significance 
threshold of 4cr and the default wavelet scales (minimum 2 
pixels, maximum 8 pixels, with a pixel size of 4"). 


3.2. Matching sources detected in more than one band 

Sources detected in more than one band were identified 
with the likelihood ratio (LR) technique. Given a distance 
r between two candidate counterparts, normalised by the 
uncertainty on the position, LR(r) is defined as the ratio 
between the probability of having a real counterpart at r 
over the probability of having a spurious counterpart at r 
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(|Pineau et al.l [2 QTi[ iFI. Similarly, the association reliability 
is defined as the probability P(i7true|^) of tho association 
being true (Ptme) conditioned on r (t he relationship be- 
tween LR and reliability is explicited in iPinean et al.ll2QTll 
eq. 11 and Appendix C). 

While the LR is usually applied to search for counter¬ 
parts in independent bands (e.g., optical counterparts of 
X-ray sources), here the bands are not independent (we ^ 
compare the 0.5-8 keV list to the 2-8 keV one; and the 0.5- ^ 

8 keV list to the 0.5-2 keV one). In this setting, and taking S 
the latter case as an example, we are testing if the position ^ 
of the source in both bands is close enough (within errors) Q 
to declare that they are the same source, against the possi¬ 
bility that the source of the additional 2-8 keV photons is 
significantly different from that of the 0.5-2 keV photons. 

In a shallow survey such as XMM-ATLAS, 

the main advantage in using the LR method over a 
nearest-neighbour match is that LR gives a likelihood for 
the match, which may help in follow-up studies, e.g. when 
searching for optical counterparts. 

The match candidates were selected on the basis of their 
positional errors; the probability of association (items [2U 
and [55] in the list in Sect. lU) was calibrated by estimat¬ 
ing a spurious LR histogram. Source IDs were assigned at 
this stage; matched sources were assigned a single ID and 
considered as a single sources from here onwards. 

3.3. Second detection pass, general settings 

In the second pass, we used the SAS EMLDetect program 
to validate the detection, refine the coordinates and obtain 
maximum-likelihood estimates of the source counts, count 
rates and fluxes. 

EMLDetect is at its core a PSF fitting code . 
Originally developed for ROSAT (|Crnd(^ce et al.l Il988l: 
iHasinger et al.l II993I ). its current versior0 has been im¬ 
proved and optimized for XMM-Newton. In particular, we 
mention the ability to operate on several individual, over¬ 
lapping pointings at the same time, and properly account 
for variations in PSF shape and in vignetting (as the tele¬ 
scope is moved between pointings, a source may be present 
in more than one pointing, but may fall on different de¬ 
tector coordinates). Rather than working on a single image 
containing the entire mosaic (like ewavelet), EMLDetect 
runs on images of individual pointings (see Sect. (Til) . 

Although EMLDetect can run on different energy bands 
at the same time, we did not use this feature, because of 
the large number of pointings and of EMLDetect limitations 
(see below), but rather we ran it three times, one per each 
band. The bands are the same as for ewavelet. In all runs, 
we used the same input list of sources, made of all matched 
sources and all unmatched sources, regardless of the bands 
in which they were detected in the first pass. 

The EMLDetect minimum likelihood was set at L = 4.6, 
as in RI3, which corresponds to a false-detection proba¬ 
bility of 1.01 X 10“^. Together with the 4cr threshold for 
ewavelet, for the final catalogue this yields a joint signifi¬ 
cance between Aa and 5cr, but which cannot be further con- 

^ iPineau et al ] (iMih also developed an Aladin plu¬ 
gin, which we used for this paper, available at 
saada.u-strasbg.fr/docs/fxp/plugin/ 

^ The XMM SAS reference manual for EMLDetect; 
xirnn.esac.esa.int/sas/13.0.0/doc/emldetect/nodeS.html 



Fig. 3. Detection grid for EMLDetect. The grey circles have 
a radius of 15^ and show the FOV of the individual pointings. 
The dashed red lines show the grid. The cyan rectangle high¬ 
lights one of the detection cells. The filled yellow circles are the 
pointings which intersect the cell, and which are used for the 
detection. Only the sources whose coordinates fall into the cell 
(the rectangle) are kept. This process is repeated for all cells 
and the final catalogue is the union of the sources from all cells. 


strained without simulations (see discussion in RI3, where a 
joint significance of 4.8cr was estimated; however this num¬ 
ber should not be immediately applied to XMM-ATLAS 
because of the different statistical properties of the back¬ 
ground, due to the different depth of this survey vs. the 
XMM-CDFS). 

The count rate to flux conversion factors were derived 
assuming a power-law spectrum with photon index F = 1.7 
and Galactic absorption of Nh = 2.3 x 10^^ cm“^, and 
by weighting the responses of the MOSI, MOS2 and PN 
cameras; they are 5.83 x 10“^^, 3.25 x 10“^^ and 1.17 x 
10“^^ erg cm“^ for the 0.5-8, 0.5-2 and 2-8 keV bands, 
respectively. 

The 8 keV upper threshold was chosen to improve the 
signal/noise ratio of detections. While quantities such as 
counts, rates and fluxes are quoted in the catalogue for the 
0.5-8 and 2-8 keV bands, fluxes in the 0.5-10 and 2-10 
keV bands can be immediately obtained by assuming the 
same F = 1.7 model spectrum, which yields the following 
conversion factors for fluxes F\ F(0.5 — I0)/F(0.5 — 8) = 
1.13 and F(2 - I0)/F(2 - 8) = 1.20. 

All EMLDetect parameters, except for what mentioned 
above (coordinate fitting, likelihood threshold, energy 
bands and conversion factors) were left at their default val¬ 
ues. The inputs to EMLDetect are images, exposure maps 
(see Sect. EH, background maps (see Sect. 13.51) and the list 
of candidate sources from Sect. 13.21 


^ Using the parameter f itposition=yes. 
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3.4. Splitting a mosaic in individual pointings 

To use EMLDetect, we had to split each mosaic-mode ob- 
sid event file into the individual pointings; for this purpose 
we used the SAS tool emosaic_prep. This tool did not set 
the RA_PNT and DEC_PNT keywords in the individual 
pointings, which were needed for further processing; how¬ 
ever, the coordinates were obtained by inspecting the atti¬ 
tude files. 

Images and exposure maps were extracted from each 
pointing’s event file in the same bands used for ewavelet. 
The SAS task eexpmap was used to compute the exposure 
maps; they are in units of s, and include corrections for 
vignetting and bad pixels. Non-vignetted exposure maps 
were also computed, to be used when producing background 
maps fSect. 1X5]) . 

Since the PSF shapes and vignetting are very sim¬ 
ilar among the three XMM-Newton EPIC cameras, for 
each pointing we summed together the images from MOSl, 
MOS2 and PN; then we did the same for the exposure and 
background maps (see below). 

3.5. Background maps 

Background maps for the individual pointings were ob¬ 
tained with__the_sa^^ used by the XMM-COSMOS 

team (|CaDDellnti et al.l I2QQ9I 1 . All input sources were ex¬ 
cised from the pointing images, producing so-called “cheese 
images”, in a number of one per pointing per energy band. A 
model made of a flat (i.e. non-vignetted, but including chip 
gaps and dead pixels) and a vignetted component (repre¬ 
senting the particle and cosmic backgrounds, respectively) 
was fit to the cheese images to obtain the background maps. 

3.6. Running EMLDetect on a wide mosaic 

EMLDetect can detect sources on a limited number of over¬ 
lapping pointings; however, the number of individual point¬ 
ings in XMM-ATLAS is too large to be analysed together. 
Therefore, we divided the ATLAS field in a grid. For each 
cell in the grid (for example, the rectangle in Fig. [3]), we 
identified the pointings whose FOV intersected the cell 
(filled circles in Fig. [3]), and ran EMLDetect on them to¬ 
gether. In this way, the number of pointings was within the 
limit allowed by EMLDetect. From the resulting detection 
list, we only kept the sources whose coordinates fell into the 
cell. We checked that no source appears twice, or was lost, 
by searching for source pairs within 3" from the cracks. We 
repeated this procedure for each of the 0.5-8, 0.5-2, and 
2-8 keV bands. The final catalogue is the union of the de¬ 
tections fro m all cells. This appr oach is similar to the one 
adopted bv iLaMassa et al.l (|2Q13[ ). 

The programme which runs EMLDetect in the cells, 
griddetect, is available (see Appendix [A]). 

3.7. Inspection of close groups 

Very close groups of sources with separation much shorter 
than the PSF width may arise in particular conditions be¬ 
cause EMLDetect is fitting the source positions and the de¬ 
tection runs are done separately for the three energy bands. 
For example, let A and B be two sources from the input cat¬ 
alogue which are present in the same EMLDetect detection 



137 136 135 

RA-TAN (deg) 


Fig. 4. Positions of detected sources, superimposed on the 0.5- 
8 keV image. 


ce 10. If at the input location of A there are not enough 
counts for a significant detection by EMLDetect, the pro¬ 
gramme may try and fit the coordinates until incorrectly 
assigning them to B’s position; this will then prevent B 
from getting a detection at the same position and in the 
same band. This might happen in one energy band but 
not in another, depending on the details of the two sources 
(e.g., input coordinates, location of the photons, hardness 
ratio). Therefore, the same EMLDetect coordinates would 
be shared by sources A and B, while their ewavelet coor¬ 
dinates would still be different. 

In order to screen and correct for this effect, we have 
looked for close groups of sources whose EMLDetect coor¬ 
dinates were closer than their 5a error. We identified 27 
groups containing 47 sources in total, and visually inspected 
all of them. For 22 groups, we shuffled the source ids to have 
the EMLDetect positions match the ewavelet ones. Five 
groups looked like genuine different sources in crowded ar¬ 
eas. 


3.8. Systematic error on astrometry 

To find whether there is any residual systematic error in 
the XMM-ATLAS astrometry, we cross-c orrelated the fi¬ 
nal c atalogue with the SDSS QSOs (DR7; ISchneider et al.l 
|2Q1Q[ ). We found a bulk shift of 0.83" in RA and —0.29" in 
DEC (XMM-ATLAS coordinate minus SDSS coordinate; 
root mean square deviation = 2.1"). We subtracted the 
above values from the XMM-ATLAS coordinates (“RA” and 
“DEC” columns only, see next Sect.), which should therefore 
not present any further shift. 


® The detection cell size is 10 x 10 pixels^ 40" x 40". 
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4. X-ray catalogue 

The XMM-ATLAS catalogue includes 1700, 1582 and 814 
sources detected by EMLDetect in the 0.5-8, 0.5-2 and 2-8 
keV bands, respectively. The number of unique sources is 
1816. The flux limits, defined as the flux of the faintest de¬ 
tected sources, are 2 x 10“^^, 6 x 10“^^ and 9 x 10“^^ erg 
s“^ cm“^ in the 0.5-2, 0.5-8 and 2-8 keV bands, respec¬ 
tively. The positions of the detected sources are shown in 
Fig. m superimposed on a 0.5-8 keV image of the mosaic. 

In addition, we list a number of 175, 103 and 47 sources 
detected in the above bands by ewavelet but not confirmed 
by EMLDetect; in the following we refer to them as sup¬ 
plementary sources. The number of unique supplementary 
sources is 234. 

All coordinates are in the J2000 reference system. The 
RA and DEC columns ar e registered to the refer ence frame 
of the SDSS-DR7 QSOs (|Schneider et al.|[2010f) . All other 
coordinate columns are not registered and present a shift; 
see Sect. 13.81 

The catalogues are available in electronic form from the 
Centre de Donnees Astronomiques de Strasbourg (CDS), 
and from the XMM-ATLAS websit^Zl, where we will publish 
updates should they become available. 

The column description is as follows; the symbols (ML) 
and (W) mark whether the column was derived from running 
EMLDetect or ewavelet, respectively. 

1. IAU_IDENTIFIER — source identifier following 
International Astronomical Union conventions; 

2. ID — unique source number; 

3-5. ID058, ID052, ID28 

6. RA — (ml) right ascension (degrees) from the 0.5-8 
keV band, if available; else, from the 0.5-2 keV band, 
if available; else, from the 2-8 keV band. This column 
has been corrected for astrometry and registered to the 
SDSS-DR7 QSO framework; 

7. DEC — (ML) declination (degrees) as above. This 
column has been corrected for astrometry and registered 
to the SDSS-DR7 QSO framework; 

8. RADEC_ERR — (ML) error on positiorjl (arcsec; 
Icr); 

9. WAV_RA — (W) merged right ascension (degrees) 
from the likelihood ratio; 

10. WAV_DEC — (W) merged declination (degrees); 

11. WAV_RADEC_ERR — (W) merged error on posi¬ 
tion (arcsec; la) 

12. WAV_RA058 — (W) right ascension (degrees) in the 
0.5-8 keV band; 

13. WAV_DEC058 — (W) declination (degrees) in the 
0.5-8 keV band; 

14. WAV_RA052 — (W) right ascension (degrees) in the 
0.5-2 keV band; 

15. WAV_DEC052 — (W) declination (degrees) in the 
0.5-2 keV band; 

^ xraygroup.astro.noa.gr/atlas.html 

^ This column was obtained by dividing the EMLDetect 
RADEC_ERR by a/2- The RADEC_ERR from EMLDetect is 

computed as (RA_ERR^-h DEC_ERR^) However, when 
two normalised one-dimensional gaussians of sigma s are com¬ 
bined, the corresponding normalised bi-dimensional gaussian 
also has sigma s, not a/2s, which is what one would get with the 
expression above. Therefore, we divided the EMLDetect value by 

V2. 


16. WAV_RA28 — (W) right ascension (degrees) in the 
2-8 keV band; 

17. WAV_DEC28 — (W) declination (degrees) in the 
2-8 keV band; 

18. RA058 — (ml) right ascension (degrees) in the 0.5-8 
keV band; 

19. DEC058 — (ML) declination (degrees) in the 0.5-8 
keV band; 

20. RA052 — (ml) right ascension (degrees) in the 0.5-2 
keV band; 

21. DEC052 — (ML) declination (degrees) in the 0.5-2 
keV band; 

22. RA28 — (ml) right ascension (degrees) in the 2-8 
keV band; 

23. DEC28 — (ml) declination (degrees) in the 2-8 keV 
band; 

24. ASSOC_RELIAB058052 — association reliabil¬ 
ity between the (WAV_RA058, WAV_DEC058) and 
(WAV_RA052, WAV_DEC052) coordinates; 

25. ASSOC_RELIAB05828 — association reliability 
between the (WAV_RA058, WAV_DEC058) and 
(WAV_RA28, WAV_DEC28) coordinates; 

26-28. SCTS058, SCTS052, SCTS28 — (ML) sum of the 
net source counts from MOS1+MOS2+PN in the 0.5- 
8, 0.5-2 and 2-8 keV bands, respectively; 

29-31. SCTS_ERR058, SCTS_ERR052, SCTS_ERR28 

— (ML) errors on SCTS058, SCTS052, SCTS28 (la); 
32-34. RATE058, RATE052, RATE28 — (ML) net count 

rates in the 0.5-8, 0.5-2, 2-8 keV bands, averaged over 
the three cameras; 

35-37. EXP_MAP058, EXP_MAP052, EXP_MAP28 — 
exposure times in the 0.5-8, 0.5-2, 2-8 keV bands, 
summed over the three cameras; 

38-40. BG_MAP058, BG_MAP052, BG_MAP28 — (ML) 
background counts/arcsec^ in the 0.5-8, 0.5-2, 2-8 keV 
bands, summed over the three cameras; 

41-43. FLUX058, FLUX052, FLUX28 — (ML) flux in the 
0.5-8, 0.5-2, 2-8 keV bands (erg s“^ cm“^); 

44-46. FLUX_ERR058, FLUX_ERR052, FLUX_ERR28 

— (ML) error on FLUX058, FLUX052, FLUX28 (la); 
47-49. DETML058, DETML052, DETML28 — (ML) detec¬ 
tion likelihoods in the 0.5-8, 0.5-2, 2-8 keV bands; 

50-51. EXT058, EXT052 - (ML) source extent in the 0.5-8 
and 0.5-2 keV band^ (a of Gaussian model in pixels; 1 
pixel = 4"); 

52-53. EXT_ERR058, EXT_ERR052 - (ML) error on 
EXT058, EXT052 (la); 

54-55. EXT_ML058, EXT_ML052 - (ML) likelihood of ex¬ 
tent in the 0.5-8 and 0.5-2 keV bands; 

56. HR — hardness ratio, computed from S = SCTS052 
and H = SCTS28 as RR = {H - S)/{H + S)- 

57. HR ERR — error on HR (la; see below); 

The columns from 58 to 87 contain source properties 
(counts, count rates, fluxes, exposure times, background, 
wavelet detection scale, source extent) from ewavelet; 
while we report them for all sources, they are actually in¬ 
teresting only for the supplementary sources. 

58-60. WAV_SCTS058, WAV_SCTS052, WAV_SCTS28 

— (W) sum of the net source counts from 
MOS1+MOS2+PN in the 0.5-8, 0.5-2 and 2-8 
keV bands, respectively; 

^ EXT28 and related columns are not included since no source 
was found to be extended in the 2-8 keV band. 


Article number, page 6 offTOl 






P. Ranalli et aL: The XMM-ATLAS survey 


61-63. WAV_SCTS_ERR058, WAV_SCTS_ERR052, 
WAV_SCTS_ERR28 — (W) errors on SCTS058, 
SCTS052, SCTS28 (la); 

64-66. WAV_RATE058, WAV_RATE052, 

WAV_RATE28 — (W) net count rates in the 0.5- 
8, 0.5-2, 2-8 keV bands, averaged over the three 
cameras; 

67-69. WAV_EXP_MAP058, WAV_EXP_MAP052, 
WAV_EXP_MAP28 — (W) exposure times in the 
0.5-8, 0.5-2, 2-8 keV bands, summed over the three 
cameras; 

70-72. WAV_BG_MAP058, WAV_BG_MAP052, 

WAV_BG_MAP28 — (W) background counts/arcsec^ 
in the 0.5-8, 0.5-2, 2-8 keV bands, summed over the 
three cameras; 

73-75. WAV_ELUX058, WAV_ELUX052, 

WAV_ELUX28 — (W) fluxes in the 0.5-8, 0.5-2, 
2-8 keV bands (erg s~^ cm“^); 

76-78. WAV_ELUX_ERR058, WAV_ELUX_ERR052, 
WAV_ELUX_ERR28 — (W) errors on 

WAV_ELUX058, WAV_ELUX052, WAV_ELUX28 

(lo-); 

79-81. WAV_WSCALE058, WAV_WSCALE052, 

WAV_WSGALE28 — (W) wavelet detection scale 
(pixels; 1 pixel = 4"); 

82-84. WAV_EXTENT058, WAV_EXTENT052, 

WAV_EXTENT28 — (W) source extent (pixels); 

85-87. WAV_EXT_ERR058, WAV_EXT_ERR052, 
WAV_EXT_ERR28 — (W) errors on 

WAV_EXTENT058, WAV_EXTENT052, 

WAV_EXTENT28 (la; pixels). 


The error on the hardness ratio (column [57l HR_ERR) 
is defined as 


HR ERR = 2 


{H + 5)2 


( 1 ) 


where H and S are the hard (SCTS28) and soft net counts 
(SGTS052), respectively, and an (SCTS_ERR28) and as 
(SCTS_ERR052) are their errors. 


5. Spectra 

To investigate the effects of the assumption of a single 
model spectrum (power-law with T = 1.7, see Sect. [3]), we 
extracted spectra for 555 sources with more than 40 counts 
in the 0.5-8 keV bancf^. 

The extraction regions are circles, whose radii were cho¬ 
sen by maximizing the signal/noise ratio for each source. 
The background regions are annuli, of inner and outer radii 
equal to 1.5 and 2 times the circle radius, respectively. In 
the case of close sources, the overlapping area was excised 
from the source and background regions. The programme 
which defines the regions, autoregions, is available (see 
Appendix 

Spectra were extracted using the cdfs-extract tool 
(see Appendix El, which takes care of identifying the 
available combinations of event files and source re¬ 
gions. cdfs-extract is a wrapper around the SAS tools 

The number of sources with SCTS058> 40 is 569, but 14 were 
dropped because they are in crowded areas and their spectra 
would include contributions from neighbour sources. 



Photon index (P) 


Fig. 5. Weighted normalised histogram of best-fit power-law 
slopes for the XMM-ATLAS sources. The inverse of the fit C- 
statistics were used as weights. 


evselect, rmf gen and arf gen to extract spectra and com¬ 
pute the response and ancillary matrices. The MOS and 
PN spectra were then summed and the response matrices 
averaged using the ETOOLs mathpha, addrmf and addarf. 

The spectra were automatically analysed with the same 
method and software (automatic XSPEC fits of unbinned 
spectra using G-sta tistic) used for the 3XMM-DR4 sources 
(|Gorral et al.ll2Ql3 ). A simple power-law model was fit to 
the data and used to derive the fluxes in the 0.5-10 and 2- 
10 keV bands. The fit was successful for 446 source^ The 
mean P is 1.7 ±0.6, consistently with the model assumed in 
Sect, m The histogram of best-fit slopes is shown in Eig. 0 
The presence of a number of flat- and inverted-spectra ob¬ 
jects can be noticed; e.g., there are 56 sources (13% of 446) 
with r < 1.0. 

The comparison between the 0.5-10 keV fluxes from the 
spectral fits and from the catalogue (the latter ones con¬ 
verted from 0.5-8 to 0.5-10 keV) is shown in Eig. [HI The 
distribution of the points around the 1-to-l line is some¬ 
what asymmetrical for EMLDetect fluxes ^ 7 x 10“^^ erg 
s“^ cm“^, with a larger dispersion above the line than be¬ 
low; this is likely due to the combination of the distribution 
of spectral slopes with the XMM-Newton effective area. 

A linear fit to the logarithms of fluxes shown in Eig. [HI 
obtained by imposing a slope of 1, yielded a normalisation: 

Espec = 1 X Feat ± 0.039 ± 0.015 (2) 

corresponding to the fluxes from spectral fits being on 
average (9 ± 4)% brighter than those from the catalogue. 
A dispersion s = 0.31 may be given as the estimate of the 
standard deviation: 

^ ~ ]\l — jy ' V (Log Fgpec ~ Log Feat) (3) 

where z/ = 1 is the number of free parameters and N is the 
number of data points, Fgpec is the flux from the spectral 
fit and Feat that from the catalogue. 

6. LogN-LogS 

The sky coverage (sky area as function of flux) for the 
XMM-ATLAS survey in the 2-8 keV band was derived with 

For the remaining 109 sources no constraint on the power-law 
slope could be put. These sources are on average fainter and 
detected with lower likelihood than those for which the fit was 
successful. 
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Flux (emldetect, erg s"^ cm“^) 



Fig. 6. Comparison of fluxes from EMLDetect with fluxes from 
spectral fits. The blue line shows the 1-to-l relationship. The 
fluxes are for the 0.5-10 keV band. 



2-10 keV flux (erg s“^ cm“^) 

Fig. 7. Coverage in the 2-10 keV band. 


Fig. 8. XMM-ATLAS Log A'-Log S in the 2-10 keV band. 
Filled circles: fluxes from the spectral fits are used whenever 
possible; catalogue fluxes otherwise. Open circles: only fluxes 
from the catalogue are used; a small horizontal shift has been 
introduced in the plot to better distinguish the error bars. The 
solid line shows the best-fit to t he filled data points . The d ashed 
line is the best-fit model from iGeorgakakis et ^ ()2008l i . The 
error bars show the la uncertainty. 

where S is the central flux of the bin; AS is the bin width; 
and for each source i, Ai is the coverage at the source’s 
flux. The sum is performed on all sources with flux Si G 
[S—AS/2, S-\-AS/2]. The error on the counts in any bin can 
be computed by assumin g ganssia n ity or , for bins with less 
than 50 sources, with the iGehrelsI (| 19861 ) approximation. 

In Fig. [8] we show the 2-10 keV differential Log A-Log S 
of the ATLAS sources detected in the 2-8 keV band. Two 
estimates are shown, from two different sets of fluxes: 

1. “spectral Log A-Log S'”: fluxes from the spectral fits 
(Sect. [5]) for 387 sources with hard detection, spectrum, 
and best-fit flux larger than 10“^^ erg s“^ cm“^(the 
threshold chosen for the faintest bin); plus fluxes from 
the catalogue (converted to the 2-10 keV band) for the 
remaining sources (i.e. those without spectra or without 
a successful fit); 

2. “catalogue Log A-Log S”: fluxes from the catalogue for 
all the 814 hard sources (also converted to the 2-10 keV 
band). 


the SAS esensmap programme, which creates a sensitivity 
map by computing count rate upper limits for each pixel. 
The exposure and background maps described in Sect. [3] 
were used as input. The coverage was then converted to 
the 2-10 keV band (Fig.[7|). 

We compute the differential counts by binning the 
sources according to their fluxes: 


n{S) 


— A — 

AS ^ Ai 

I 


(4) 


It can be noticed that the spectral Log A-Log S extends 
to brighter fluxes than the catalogue one, by a factor < 2, 
which is consistent with the use for a few sources of a F 
flatter than 1.7. 

F or comparison, we a l so plo t in Fig.[8]the best-fit model 
from IGeorgakakis et all (|2008[ ) , which was derived from a 
joint analysis of several surveys, some of which are much 
deeper than XMM-ATLAS. The model has a broken power- 
law shape, with a break flux at 1.2 x 10“^^ erg s“^ cm“^, 
very close the XMM-ATLAS flux limit (at 1.3 x 10“^^ erg 
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cm“^ the coverage is 5.3% of the nominal area, and be¬ 
comes ill-defined at fainter fluxes). Consistently, the XMM- 
ATLAS Log N-Log S does not show such a break. 

A linear fit (weighted least squares) to the spectral 
Log A^-Log S yielded LogA^ = (—2.47 ± 0.05)Log5' — 
18.0 T 0.7 (errors at la ), consistent within errors with 
iGeorgakakis et all (|2008[ ) and with a Euclidean distribu¬ 
tion. A fit to the catalogue Log N-Log S yielded LogN = 
(—2.37 ±0.08)LogS' —20.09 ± 1.03, which is only consistent 
with Euclidean counts at a 2cr level. 


7. Conclusions 

We have presented the observations, data reduction, cat¬ 
alogue and number counts of the XMM-ATLAS survey, 
which covers 7.1 deg^ with flux limits (defined as the flux 
of the faintest detected sources) of 2 x 10“^^, 6 x 10“^^ and 
9 X 10“^^ erg s“^ cm“^ in the 0.5-2, 0.5-8 and 2-8 keV 
bands, respectively. 

We derived the catalogues with a two-step procedure. 
Eirst, the ewavelet SAS task was used to identify candidate 
sources with a significance equivalent to 4cr, and to find 
their coordinates. Next, we used the SAS EMLDetect tool 
to further check the significance of the sources, and obtain 
counts, count rates and fluxes. The final catalogues contain 
1700, 1582 and 814 sources in the 0.5-2, 0.5-8 and 2-8 keV 
bands, respectively, with a total of 1816 unique sources. A 
list of supplementary sources, detected by ewavelet but 
not confirmed by EMLDetect, is also provided. 

To investigate the effect of assuming a common spectral 
model for all sources to convert count rates to fluxes, we 
extracted spectra for each source with at least 40 counts 
and fitted them with a power-law model. We found that 
on average, the fluxes from the spectral fits are (9 ± 4)% 
brighter than those from assuming a common power-law 
photon index of T = 1.7. Also, the average best-fit spectral 
slope is r = 1.7 ± 0.6. 

We derived the 2-10 keV differential Log N-Log S for 
the XMM-ATLAS sources, which spans the 10“^^-10“^^ 
erg s“^ cm“^ flux interval. Using the fluxes from the spec¬ 
tral fits (for sources with spectra; and fluxes from the cata¬ 
logue for all other sources) produces a Log A^-Log S which 
is less noisy and more consistent with Euclidean counts 
than using fluxes from the catalogue. A weighted linear fit 
yielded LogA^ = (—2.47 ± 0.05)LogS'— 18.0 ± 0.7 (errors at 
la); this is consistent within errors with previous studies, 
e.g. iGeorgakakis et al.l (|2008h . 

Einally, we release the software tools which have been 
developed or enhanced to accomplish the analysis in 
this paper: griddetect, a programme to run EMLDetect 
over a very wide mosaic; autoregions, to define extrac¬ 
tion regions for spectra and aperture photometry; and 
cdfs-extract, to extract spectra for multiple sources from 
multiple observations. 

Appendix A: Released software 

Eor the software described below, more information and 
instructions for use can be found on the author’s websit43 
and source repositorAl^. 

members.noa.gr/piero.ranalli 

github.com/piero-ranalli 


cdfs-extract and autoregions 

Given a number of event files and a list of source and back¬ 
ground positions, the cdf s-extract programme checks if 
the source and background are in the field of view of any 
observation, and it extracts products accordingly. The ob¬ 
servations may or may not overlap. 

The extracted products are spectra (source and back¬ 
ground) and responses (source and, optionally, background 
RMEs and AREs), or aperture photometry, or lightcurves. 
The programme takes as input two text files, containing: 

— a list of event files, exposure maps and (optional) im¬ 
ages; 

— a list of source and background positions and radii (the 
camera may optionally be specified, allowing the same 
source to have different positions / radii for different cam¬ 
eras) . 

Spectra and responses can be summed if the user wants, 
producing a single spectral file for each source and each 
camera. 

The list of source and background positions can be ei¬ 
ther be manually written, or automatically produced by 
autoregions. If the sources are not too close to each other 
and the background has no large variations on scales < 20- 
30" then annuli can be used as background regions and their 
positions and radii can be automatically generated (i.e., the 
user is dealing with a shallow survey such as XMM-ATLAS, 
XXL or XMM-COSMOS; conversely, deep surveys such as 
the XMM-CDES feature close groups of sources and com¬ 
plex spatial patterns in the background which make auto¬ 
matic procedures unreliable). 

Given a catalogue of sources, stored in a database 
(e.g.: PostgreSQL) and containing source counts and back¬ 
ground surface brightnesses (i.e., the EMLDetect output), 
the autoregions programme computes extraction regions 
by maximizing the expected signal/noise ratio. Overlapping 
regions are identified, and where possible, the overlapping 
areas are excised. 

Besides this work, the programmes described above have 
been, or are being used, in several papers, e.g. t he obscured 
AGN in XMM-COSMOS (iT.anznisi et, the XMM- 

CDES spectral survey (Comastri et ah, in prep.), the XXL 
brightest AGN spectral survey (Eotopoulou et ah, in prep.). 

Both cdf s-extract and autoregions are written in 
(modern-style) Perl using the Perl Data Language (PDL) 
libraries, and are free software released under the terms of 
the GNU Affero GPL licens43- 

griddetect and associated libraries 

The motivation for griddetect, and how it works, have 
already been described in Sect. 13.61 (see also Eig. [3]). 

Here we also mention that the XMM-Newton SAS pro¬ 
vides a tool (emosaicproc), whose functionality is partly 
duplicated by griddetect. The main difference between 
the tools may be summarized as follows: emosaicproc au¬ 
tomatically computes the grid, and requires less data prepa¬ 
ration by the user if a single mosaic-mode obsid is to be 
processed, but only works on a single obsid; griddetect 
needs the user to define the grid cracks, but it is easier to 

The Affero GPL is a variant of the GPL which also adds 
protection when the software is used as a web service. 
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add pointings from different obsids, and it can also compute 
background maps. 

griddetect works on individual pointings, splitted from 
mosaic-mode ev ent files; it also compute s background maps 
according to the lCappellnti et al.l (|2QQ9[ ) method. The input 
comprises: 

— the grid specification: a rotation angle, and the x and y 

coordinates (rotated RA and DEC) of the grid lines; 

— the list of pointings with their coordinates. 

The output is a series of two catalogues per detection cell: 
one containing the detected sources, and another one con¬ 
taining only the sources within 3" from the cracks (used to 
check if there are missing or repeated sources, arising be¬ 
cause of numeric effects placing them alternatively on the 
two sides of the crack). 

Together with griddetect, we also release a set of as¬ 
sociated libraries which may be of more general purpose: 
XMMSAS: :Extract and XMMSAS: :Detect. These are object- 
oriented Perl packages which present a high-level interface 
to some SAS and FTOOLS commands. They can be called 
by Perl programmes to filter event files and extract images 
and exposure maps using either pre-defined filters (e.g., the 
instrumental lines excluded in Sect. HI) or user-defined ones; 
and to compute background maps, source masks and call 
EMLDetect. 

griddetect and its associated libraries are written in 
(modern-style) Perl using the Moose object system and the 
Perl Data Language (PDL) libraries, and are free software 
released under the terms of the GNU Affero GPL license. 
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